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Abstract. Recent cosmological observations hint at a deviation from the simple power-law 
form of the primordial spectrum of curvature perturbations. In this paper we show that in 
the presence of a tensor component, a turn-over in the initial spectrum is preferred by current 
observations, and hence non-power-law models ought to be considered. For instance, for a 
power-law parameterisation with both a tensor component and running parameter, current 
data show a preference for a negative running at more than 2.5er C.L. As a consequence 
of this deviation from a power-law, constraints on the tensor-to-scalar ratio r are slightly 
broader. We also present constraints on the inflationary parameters for a model-independent 
reconstruction and the Lasenby <fe Doran (LD) model. In particular, the constraints on the 
tensor-to-scalar ratio from the LD model are: tld = 0.11 ± 0.024. In addition to current 
data, we show expected constraints from Planck-like and CMB-Pol sensitivity experiments 
by using Markov-Chain-Monte-Carlo sampling chains. For all the models, we have included 
the Bayesian Evidence to perform a model selection analysis. The Bayes factor, using current 
observations, shows a strong preference for the LD model over the standard power-law pa- 
rameterisation, and provides an insight into the accuracy of differentiating models through 
future surveys. 
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1 Introduction 

Inflationary models have the merit that they not only explain the homogeneity of the universe 
on large-scales, but also provide a theory to explain the origin of perturbations as observed in 
the Cosmic Microwave Background (CMB). During the inflationary period, quantum fluctu- 
ations of the field were driven to scales much larger than the Hubble horizon and eventually 
turned into density perturbations (scalar) observed in the CMB and the large-scale dis- 
tribution of galaxies, together with a gravitational wave (tensor) contribution. The vector 
contributions, however, are expected to be negligible since these modes decayed very rapidly 
once they entered the Hubble horizon. The scalar and tensor contributions are summarised 
by the primordial power spectrum Viz(k) and Vj-(k), respectively. The scalar spectrum for 
a single-inflaton field <fi, in the slow-roll approximation, is given by [21]: 



where the expression is evaluated at the horizon exit k = aH. Since the initial spectrum is an 
unknown function, one needs to carry out a full numerical calculation from the onset of the 
inflationary phase, or to assume a particular functional form of it. The simplest proposal is 
to parameterise the shape of T > Ti(k) by a power-law. Although the power-law assumption has 
provided reasonable agreement with cosmological observations, some recent analyses have 
shown that if a running of the scalar spectral- index is taken into account, there exists a pref- 
erence for a negative running-value at 1.8<7 C.L. from WMAP7+ACT (SPT) measurements 
[10, 15] and at 2.2cj C.L. with WMAP7+QUaD [4]. It has also been shown that the existence 
of a turn-over in Vn(k), by using model-independent techniques, is preferred [13, 14, 23, 31]. 
The presence of this turn-over plays an important role in explaining current cosmological 
observations and cannot therefore be ignored when constraining the inflationary parameters. 
For instance, in the slow-roll approximation, the shape of the spectrum of tensor perturba- 
tions is 




(1.1) 



V T ik)= -H 2 
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which depends on the form of the scalar spectrum, and vice- versa, via the potential of the sin- 
gle scalar-field. To place constraints on the amplitude of tensor contributions, it is customary 
to define the tensor-to-scalar ratio as 



The dependence of the scalar spectrum on the tensor spectrum is evident in the Lasenby & 
Doran model, where both spectra depend upon the same best-fit parameters. In a previous 
paper [31], we found that standard ACDM models with a turn-over in the scalar spectrum are 
preferred over a simple power-law parameterisation. In this work, by assuming a power-law 
parameterisation of the tensor spectrum, we show that the bending of the scalar spectrum 
is enhanced due to the presence of a tensor component. To avoid misleading results due to 
the particular choice of parameterisation, the shape of the scalar spectrum is described by 
employing a model-independent reconstruction. We then show that current constraints on 
the tensor-to-scalar ratio (1.3) are broadened for non-power law V-R,(k) models. We also dis- 
cuss the constraints on r for a massive scalar-field in the Lasenby & Doran model. Finally, 
by considering future experiments we present their expected constraints on the inflationary 
parameters. For all the models, the Bayes factor is computed in order to perform a model 
comparison. 

The paper is organised as follows: in the next Section we list the data sets and the 
cosmological parameters considered. In Section 3 we study different models suggested to de- 
scribe the form of the scalar spectrum. Then, we show the resulting parameter constraints 
on the tensor-to-scalar ratio and the preferred form of the power spectrum using current cos- 
mological observations. In the same section we provide future constraints on r expected by 
Planck-like and CMB-Pol experiments. Performance assumptions for Planck and CMB-Pol 
are taken from [25] and [2]. We present our conclusions in Section 5. 



2 Theoretical Framework 

Even though the primary parameters in the standard ACDM model have already been tightly 
constrained and have little impact on the B-mode spectrum, it is worthwhile to perform a 
full parameter-space exploration to determine the tensor-to-scalar ratio constraints in each 
model. We assume purely Gaussian adiabatic scalar and tensor contributions in a flat ACDM 
model 1 specified by the standard parameters: the physical baryon Q\,h 2 and cold dark matter 
density Q c h 2 relative to the critical density (h is the dimensionless Hubble parameter such 
that Hq = lOO/i kms~ 1 Mpc~ 1 ), 6 is lOOx the ratio of the sound horizon to angular diameter 
distance at last scattering surface, r denotes the optical depth at reionisation. We consider 
the tensor-to-scalar ratio for each model i, as = VT(i)(,k)/V-R,(i)(k)- Hereafter, we set the 
tensor-to-scalar ratio rj = rj(A^), at a scale of k^ = 0.015Mpc _1 . A study of the appropriate 
scale to use is given by [8]. Aside from the Sunyaev-Zel'dovich (SZ) amplitude A$z used 
by WMAP analyses, the ACT likelihood incorporates two additional nuisance parameters: 
the total Poisson power A p and the amplitude of the clustered power A c . The parameters 
describing the primordial spectra for each model are listed in the next section, together with 

1 Except for the LD model, which is based on a marginally closed universe Qk < 0. 
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Figure 1: Comparison of theoretical prediction of r = 0.1 inflation model, and the WMAP, 
BICEP and QUaD data for the B-mode power spectrum. 



the flat priors imposed in our Bayesian analysis. 

Throughout the analysis, the theoretical temperature and polarisation C/s spectra are 
generated with a modified version of the CAMB code [20] , and the parameter estimation is 
performed using the CosmoMC program [19]. The calculation of the Bayesian evidence Z, to 
perform the model selection, requires a multidimensional integration over the likelihood and 
prior. To do this, we make use of the MultiNest algorithm [11, 12]. The Bayes factor Bij, 
or equivalently the difference in log evidences In Z{ — In Zj , provides a measure of how well 
model i fits the data compared to model j. In order to make a qualitative model comparison, 
we consider the Jeffreys guideline: if B{j < 1 model i should not be favoured over model j, 
1 < < 2.5 constitutes significant evidence, 2.5 < By < 5 is strong evidence, while B^ > 5 
would be considered decisive [30, 32]. 



Current Observations 

To compute posterior probabilities for each model we use temperature and polarisation mea- 
surements from the Wilkinson Microwave Anisotropy Probe 7- year (WMAP7; [16]) and the 
Atacama Cosmology Telescope (ACT; [10]) data. To improve polarisation constraints, we 
include observations from QuaD [4], whose primary aim is high resolution measurements 
(154 < I < 2026) of the E-mode signal, and BICEP data [7] which probes intermediate 
scales (21 < £ < 335). Figure 1 shows the .B-mode spectrum predicted from a power-law 
parameterisation, with r = 0.1, along with la constraints obtained by using current obser- 
vations. In addition to CMB data, and to strengthen the constraining power, we incorporate 
distance measurements from the Supernova Cosmology Project Union 2 compilation (SCP; 
[1]) and Large Scale Structure data from the Sloan Digital Sky Survey (SDSS) Data Release 
7 (DR7) Luminous Red Galaxy (LRG) power spectrum [28]. We also also consider baryon 
density information from Big Bang Nucleosyntesis (BBN, [6]) and impose a Gaussian prior 
on the Hubble parameter today Hq from measurements of the Hubble Space Telescope (HST; 
[29]) key project. 
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Figure 2: Polarization noise power spectra for forthcoming experiments. Note that these 
curves include uncertainties associated with the instrumental beam. The red line shows the 
B-mode power spectrum for the standard inflationary model with r = 0.1. 



Future surveys 

In order to forecast expected constraints of future CMB observations, one can make use 
of the Fisher information matrix under the assumption that each parameter is Gaussian- 
distributed. However to obtain constraints more closely related to what is expected by future 
experiments, we perform a Monte Carlo analysis sampling over all variables involved in the 
description of the CMB spectrum. We simulate future experiments by generating mock data 
of the Ci 's from a xt^+i distribution with variances [26]: 



(2£ + l)f sky 
2 

(2£ + l)f sky 



{c xx + N xx ) \ 



{Cj E f + (CT + Nf) (C EE + N EE ) 



(2.1) 
(2.2) 



where X = T,E and B label the temperature and polarisations; f sky is the fraction of the 
observed sky. The C^ Yl s represent the theoretical spectra and Y the instrumental noise 
spectra for each experiment. In experiments with multiple frequency channels c, the noise 
spectrum is approximated [3] by 



N, 



x 



(2.3) 



where the noise spectrum of an individual frequency channel, assuming a Gaussian beam, is 



t(£ + l) 



9 2 

fwhm 



81n2 



Sxy- 



(2.4) 



(cTpix #fwhm) 2 exp 

The pixel noise from temperature and polarisation maps are considered as uncorrelated. The 
noise per pixel a^ ix (and cr^ ix = V2(Tp ix ) depends on the instrumental parameters; #f w h m is 
the full width at half maximum (FHWM) of the Gaussian beam. 



- 4 - 



For the Planck experiment, we include three channels with frequencies (100 GHz, 143 GHz, 
217 GHz) and noise levels per beam (fJ^ ix ) 2 = (46.25 //K 2 , 36 /iK 2 , 171 /jK 2 ). The FHWM 
of the three channels are #f w hm =(9.5,7.1,5.0) arc-minute. These figures are taken from 
the values given in [25]. We combine three channels for the CMBPol experiment [2] with 
frequencies (100 GHz, 150 GHz, 220 GHz), noise levels (crj ix ) 2 = (729 nK 2 , 676 nK 2 , 1600 
nK 2 ) and 6>f w h m = (8,5,3.5) arc-minute. Sky coverages of / s k y = 0.65,0.8 are respectively 
assumed and integration time of 14 months. In Figure 2, we show the noise levels for these 
experiments as a function of multipole number I. The blue line corresponds to the B-mode 
power spectrum using the standard power- law parameterisation, with r = 0.1. The lensed 
Cf is also shown in the same Figure, which can be treated as a part of the total noise power 
spectrum NP as well as the instrumental noise power spectra [24]. For more information of 
the noise and beam profile of each frequency channel, please refer to [22]. 

3 Primordial power spectra constraints 
3.1 Power-law parameterisation 

Because slow-roll inflation predicts the spectrum of curvature perturbations to be close to 
scale-invariant, the simplest proposal is to assume that the initial spectrum has a power-law 
form, parameterised by 



where the spectral index n s is expected to be close to unity. A spectrum where the typical 
amplitude of perturbations is identical on all length scales is known as Harrison-Zel'dovich 
spectrum (n s = 1), and it has been ruled out by several studies (see for instance [31]). Here, 
we assume, for simplicity, that the tensor spectrum is also described by a power-law function: 



where the tensor amplitude At is related to tensor-to-scalar ratio r s = At/A s . For this pa- 
rameterisation, we assume that r(ko) and the tensor spectral index rit satisfy the consistency 
relation for a single field slow-roll inflation n t = — r s /8 [9]. The power-law parameterisation 
thus contains only three free parameters: A3, n s , and r s . For these parameters, we assume 
a prior A s = [1,50] x 10~ 10 for the amplitude, a conservative prior for the spectral index 
n s = [0.7, 1.2] and a tensor-to-scalar ratio prior of r s = [0, 1]. 

Figure 3 shows ID and 2D marginalised posterior distributions of the scalar spectrum 
index n s and the tensor-to-scalar ratio r s , using both current cosmological observations 
(black line) and future experiments (red for Planck and green for CMBPol). The bottom 
panel shows the limits imposed by current and future experiments. For present observations: 
n s = 0.964±0.011 and r s < 0.171 (mean values of 68% C.L. are quoted for two-tailed distribu- 
tions, whilst one-tailed distribution only the upper 95% C.L.). These results are in agreement 
with previous studies, i.e. [10, 15, 16]. With regards to future constraints, we have used mean 
values obtained from current observations as the fiducial model (with fixed r s = 0.1). We 
notice that la error bars of the spectral index n s , shown in the bottom panel of the same Fig- 
ure, reduce by about four times using a Planck-like experiment and five times for a CMBPol 
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Figure 3: Left panel: ID and 2D probability posterior distributions for the power spectrum 
parameters, assuming a simple tilt parameterisation (n s ); using both current observations 
(black line) and future experiments (red for Planck and green for CMBPol). 2D constraints 
are plotted with la and 2a confidence contours. Right panel: Reconstruction of the scalar 
spectrum using present data; lighter regions represent an improved fit. 



experiment. Whereas Planck will be able to distinguish tensor components with an accuracy 
of a r = 0.013, this is highly improved by CMBPol data a r = 0.0009. If we consider only 
one channel for comparison, e.g. 100 GHz, the constraints on the tensor-to-scalar ratio are 
given by a r = 0.02, in agreement with previous results [5]. The top-right panel of Figure 3 
illustrates the resulting shape of Vn(k) corresponding to the posterior distributions using 
present data. 



3.2 Running scalar spectral-index 

A further extension is possible by allowing the scalar spectral index to vary as a function of 
scale, such that n s (k). This can be achieved by including a second order term in the expansion 
of the power spectrum 

Vn(k) = A s (-) , (3.3) 

where n run = dn/dlnk is termed the running of the tilt and we would expect n mn ~ for 
standard inflationary models. We have kept the same tensor spectrum as in the simple power- 
law parameterisation, with a tensor-to-scalar ratio r mn at a scale of k^ = 0.015 Mpc -1 to 
avoid correlations amongst parameters [8] . We maintained the same priors for the inflationary 
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Figure 4: Left panel: ID and 2D probability posterior distributions for the inflationary pa- 
rameters, assuming a power-law with a running parameter (n run ); using both current data 
(black line) and future experiments (red for Planck and green for CMBPol). 2D constraints 
are plotted with la and 2a confidence contours. Right panel: Reconstruction of the scalar 
spectrum using present data; lighter regions represent an improved fit. The top label de- 
notes the Bayes factor of the n run -model compared to the power-law ra s -model, using current 
observations. 

parameters A3, n s , and r run and select a prior of the running parameter of n mn = [—0.1, 0.1] 
as used by [23]. 

Figure 4 shows the ID and 2D marginalised posterior distributions for the inflationary 
parameters, using current experiments (black line): n s = 0.985 ±0.017, n run = —0.043 ±0.018 
and r mn < 0.324; and Planck (red line) and CMBPol (green line) realisations. The top label 
of the figure indicates the Bayes factor using present observations, which in this case and 
throughout the paper is compared with respect to the power-law parameterisation. We first 
note that in the presence of a tensor component, the bending of the scalar spectrum is en- 
hanced through a larger running parameter 2 . We also observe that using current experiments 
a negative n run parameter in preferred by more than 2.5<r C.L. Hence the necessity to include 
a turn-over in the power spectrum. This result is confirmed by noticing the Bayes factor 
is significantly favoured compared to the simple power-law model, B niun ,n s = ±2.0 ± 0.3. 
Considerations of the running of running of the spectral index are also being explored [27]. 

2 The constraints of the running parameter without tensor components are n run = —0.028 ± 0.014. 
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We notice that correlations created by the inclusion of the running parameter broaden the 
constraints on the tensor-to-scalar ratio by about 1.5 times. Future constraints are also broad- 
ened compared to the power-law parameterisation. The summary of the constraints on the 
inflationary parameters is shown in the bottom panel of Figure 4, and the reconstruction of 
Vn(k), using present data, in the top-right panel. 

3.3 Model independent reconstruction 

We have seen that deviations from the simple power-law, by the introduction of the running 
parameter, are relevant in explaining present data. In order to corroborate this result and 
look for deviations from the power-law parameterisation, we consider a model-independent 
reconstruction. The reconstruction process we follow is based on the approach used previ- 
ously by [31]. We place two fixed A:-nodes at sufficiently separated positions [/c m i n , A: max ], with 
varying amplitudes L4 Sj fc min , ^4 s ,fc max ]) an d place inside additional 'nodes' with the freedom to 
move around in both position ki and amplitude ^4 s ,fc;- We assume that most of the astro- 
physical information is encompassed within the scales /c m ; n = 0.0001 Mpc -1 and k max = 0.3 
Mpc -1 . Outside of these limits we consider the spectrum to be constant with values equal to 
those at k m [ Q and /c max respectively. We allow variations in the amplitudes with a conservative 
prior A s ^. 6 [1,50] x 10 -10 . To maintain continuity between /c-nodes, a linear interpolation 
is performed such that the form of the power spectrum is described by 

P-nik) = < A s>ki k min < h < k i+ i < /c max , (3.4) 

. ^s,fc max k > &max> 

and with linear interpolation for k m i n < ki < k max . 

We have restricted the model-independent reconstruction to two internal-nodes which we 
consider are sufficient to provide an accurate description of the shape of the power spec- 
trum. The tensor spectrum is parameterised by a power-law form, similarly to the one in 
Section 3.2. Here, the tensor-to-scalar ratio, given by r2 ki = At/Pn(ko), is computed at the 
scale ko = 0.015 Mpc -1 , and also satisfies the consistency relation = — n t/&] with prior 
^ = [0,1]. 



The top-left panel of Figure 5 presents the ID and 2D marginalised posterior distribu- 
tions for the parameters used in the model-independent reconstruction. At the largest scales, 
we observe the lack of tight constraints on the amplitude A±, mainly due to the cosmic 
variance and correlations with other parameters. At smaller scales, the constraints on the 
amplitues (i.e. A2, A3 and A4) get tighter. We notice the presence of a bi-modal distribution 
in the medium/small scales, represented by fai-, where the highest peak (k ~ O.OlMpc -1 ) 
matches the position of the turn-over in the primordial spectrum, as seen in the top-right 
panel of Figure 5. The other peak is located where the constraints seem to improve by up- 
dated data sets: at the overlapping of WMAP/ACT observations (0.1 < k < 0.14) with 
LRG7 measurements. The reconstructed spectrum clearly presents a turn-over, however with 
the bending at small scales less pronounced that in the running model. Notice that the Bayes 
factor, shown in the top label of the same Figure, is significantly preferred over the simple 
tilt model, even though the two-internal-node reconstruction contains four additional param- 
eters; it is also marginally preferred over the running model. Future experiments will be able 
to pin-down accurately the shape of the primordial spectrum at medium and small scales 
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Figure 5: Left panel: ID and 2D probability posterior distributions for the power spectrum 
parameters, assuming a two internal-node reconstruction (2fcj); using both current cosmolog- 
ical observations (black line) and future experiments (red for Planck and green for CMBPol). 
2D constraints are plotted with la and 2a confidence contours. Right panel: Reconstruction 
of the scalar spectrum using present data; lighter regions represents an improved fit. Top 
label denotes the Bayes factor of the 2fcj-model compared to the power-law n s -model, using 
current observations. 



(^2), however at the largest scales (k\) the cosmic variance still dominates, as seen in the ID 
posterior distribution of A\. Current and future constraints of the inflationary parameters 
are summarised in the bottom panel of Figure 5. 

3.4 Lasenby & Doran model 

The Lasenby & Doran model is based on the restriction of the total conformal time available in 
a closed universe [18]. At the largest scales, the predicted scalar and tensor spectra naturally 
incorporate a drop-off without the need to parameterise them, whilst, at small scales they 
mimic a slight running behaviour. An important point to bear in mind, is that in the LD model 
the functional forms of H and <f> during inflation are expressed using just two parameters 
bo and 64 [31, 32]. These parameters describe the initial conditions, along with the standard 
cosmological parameters, and therefore the primordial spectra generated by the LD model 
are given in terms of 

Pn(k) = Pn{k; b , 64, fi., H ), P T {k) = P T (k; b , 64, SU, H ). (3.5) 
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Figure 6: Derived primordial power spectra, Pji(k) and Pj-(k), from the Lasenby &: Doran 
model, using a different set of parameters of bo and 64 (left and middle panels); units of 64 
are given in [10 -10 ]. The right panel shows the tensor-to-scalar ratio tld = Pt/Ph- 

Notice that the tensor-to-scalar ratio tld is a derived quantity in terms of the cosmological 
parameters, Ho, fli, and the initial-conditions parameters 60 and 64: 

riD(k) = r LD (k; b , h, H ). (3.6) 

That is, if we use values of bo and 64 along with the cosmological parameters there is no need 
to introduce additional variables to describe the tensor-to-scalar ratio t^d. Figure 6 shows the 
primordial spectra, both scalar and tensor, for a given combination of 60 and 64 parameters. In 
the right panel of this Figure, we illustrate the tensor-to-scalar ratio and its degeneracy with 
a selection of the parameters, for instance, the combination of {bo = 3.2, 64 = —10 x 1CP 10 } 
or {bo = 3.0,64 = —12 x 10~ 10 } provides the same tensor-to-scalar ratio, even though their 
scalar and tensor spectra differ considerably. For further details about the LD model see, for 
instance [17, 18, 32]. To compute the LD spectra we refer to [31, 32]. We have also chosen the 
priors based on the same paper: = [—0.05, 10 -4 ], bo = [1,4], 64 = [—30, —1] x 10~ 9 . Figure 
7 shows ID and 2D marginalised posterior distributions of the parameters involved in the 
description of the LD model. A novel result from the LD model is that its constraints on the 
tensor-to-scalar ratio are different from zero: tld = 0.11 ± 0.024, contrary to the standard 
power-law parameterisation. This happens mainly due to the 2 -type potential assumed in 
the model. The Bayes factor of this model compared to the simple-tilt model is shown in the 
top label of the same Figure. The low number of parameters and the reduced power at both 
large and small scales make the LD model strongly favoured compared to the simple tilt and 
significantly so compared to the running and the two-internal-node reconstruction. Future 
experiments will provide an insight on discriminating amongst models, as we will see in the 
next section. 

4 Model Selection 

Throughout the analysis, we have included the Bayes factor for each model and found that the 
Lasenby & Doran model is the most preferred by current observations. Future experiments 
will be able to distinguish between models more effectively. Let us assume for a moment that 
the LD spectra represent the true model. We then use the LD spectra, with best-fit values 
obtained by using present data (shown in the bottom panel of Figure 7), as the fiducial model 
to simulate future CMB observations. We analyse this mock data to reconstruct the input 
spectrum using the set of models aforementioned. Table 1 shows the Bayes factor for the 
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Figure 7: ID and 2D probability posterior distributions for the power spectrum parameters, 
assuming a Lasenby & Doran model (LD); using both current cosmological observations 
(black line) and future experiments (red for Planck and green for CMBPol). 2D constraints 
are plotted with la and 2a confidence contours. The top label denotes the associated Bayes 
factor with respect to the power-law n s model using present data. *In this model, r^D is a 
derived parameter. 



different models compared to the LD model. There is indeed a distinction between models, 
with the data clearly indicating a preference for the LD model, used to generate the input- 
simulated data. Idealised Planck results might provide decisive conclusions on distinguishing 
the LD model from the simple-tilt ra s -model, B nst LD = —6.3 ± 0.3, and a running n run - 
model, B nrunt LD = —6.5 + 0.3, and strong preference when compared to the two-internal-node 
reconstruction 2/cj-model, B2k it LD = —3.1 ± 0.3. There will also be a strong preference for 
the model independent reconstruction over both the n s and n run models: Z^fc^ns = +3.2 ±0.3 
and 02fcj,n run = +3.4 + 0.3, respectively. With regards to the CMBPol experiment, this might 
definitely differentiate the LD spectrum from the rest of the spectra. In contrast to the Planck 
experiment, the model-independent reconstruction for CMBPol is now strongly favoured 
compared to the n run model. CMBPol also provides a strong preference to differentiate the 
simple-tilt model n s over the running model n run : £>n s ,nrun = +2.5 ± 0.3. Therefore, future 
experiments certainly will be able to differentiate between these models and pin down the 
right form of the primordial spectrum. 
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Table 1: Model Selection. The input spectrum, given by the LD model, is reconstructed 
using different models. We show the Bayes factor for each model i, Bi t LD, compared to the 
LD model. 



Planck CMBPol 

"LD 0.0 ±0.3 0.0 ±0.3 

n s -6.3 ±0.3 -13.0 ±0.3 

n run -6.5 ±0.3 -15.5 ±0.3 

2h -3.1 ±0.3 -10.2 ±0.3 



5 Discussion and Conclusions 

In this paper we have performed a MCMC exploration of the full cosmological parameter- 
space and showed current and future constraints on the inflationary parameters, with par- 
ticular attention to the tensor-to-scalar ratio. We have considered models that deviate from 
the standard power-law in the scalar power-spectrum: a power-law parameterisation with 
running behaviour and the spectra predicted from the Lasenby & Doran model. By imple- 
menting a model- independent reconstruction for Vn(k), we found that a turn-over in the 
scalar spectrum is preferred to explain cosmological observations. The inclusion of a ten- 
sor component enhances this turn -over, for instance, Ti Tun — -0.043 ± 0.018 compared to 
n mn = —0.028 ± 0.014 without tensors. A similar form of the scalar spectrum has been pre- 
viously obtained assuming different model-independent reconstructions, and some of them 
with different data sets [13, 14, 31]. We have assumed that a power-law parameterisation of 
Vf(k) is sufficient to describe current data. Even though we have not given the results for 
the standard cosmological parameters Q^h 2 , il c /i 2 , 8, r, their best-fit values remained essen- 
tially unaffected throughout the models. For all the models, we have computed the Bayes 
factor and compared each to the simple power-law parameterisation. We found, using cur- 
rent observations, that the preferred model is given by the LD model. The summary of the 
analysis, illustrated in Figure 8, displays how the constraints on the tensor-to-scalar ratio are 
broadened for non-power law models. We notice that the best-fit value of r run is slighly offset 
from zero and coincides with the peak of tld- Throughout the models, the tensor-to-scalar 
ratio has been computed at a particular scale k$ = 0.015 Mpc -1 . However, to illustrate the 
robustness of the model selection, over a different choice of scale ko, we compute the Bayesian 
evidence for all models at k$ = 0.002 Mpc -1 . The results are esentially unaffected and still 
show a preference for the LD model. The Bayesian evidece, for each model, compared to the 
power-law parameterisation (n s ) are as follow: 

£> mn ,n s £>2ki,n s ^LD,n s 

±1.8 ±0.3 ±2.29 ±0.3 ± 3.0 ± 0.3 

With regards to future surveys, the Planck satellite will be able to differentiate the 
running and tilt model from the LD model, but not decisively from the two-internal-node 
reconstruction. The improvement using CMBPol selects the right form of the primordial 
spectrum, as shown in Table 1. 
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Figure 8: 1-D marginalised posterior distributions of the tensor-to-scalar ratio for the dif- 
ferent models (left panel), along with their Bayesian evidence and number of parameters of 
each model (right panel). The Bayes factor is compared to the simple-tilt model (n s .) 



Acknowledgments 

This work was carried out largely on the Cambridge High Performance Computing cluster, 
DARWIN. JAV is supported by CONACYT Mexico. 



References 

[1] R. Amanullah and et. al. Spectra and Hubble Space Telescope Light Curves of Six Type la 
Supernovae at 0.511 < z < 1.12 and the Union2 Compilation. The Astrophysical Journal, 
716(1):712, 2010. 

[2] B-Pol Collaboration. Exper. Astron, 23(5), 2009. 

[3] M. Bowden and et. al. Scientific optimization of a ground-based CMB polarization experiment. 
Monthly Notices of the Royal Astronomical Society, 349(l):321-335, 2004. 

[4] M. L. Brown and et. al. Improved Measurements of the Temperature and Polarization of the 
Cosmic Microwave Background from QUaD. The Astrophysical Journal, 705(1):978, 2009. 

[5] C. Burigana, C. Destri, H. J. de Vega, A. Gruppuso, N. Mandolesi, P. Natoli, and N. G. 
Sanchez. Forecast for the Planck Precision on the Tensor-to-Scalar Ratio and Other 
Cosmological Parameters. The Astrophysical Journal, 724(1):588, 2010. 

[6] S. Buries, K. M. Nollett, and M. S. Turner. Big Bang Nucleosynthesis Predictions for Precision 
Cosmology. The Astrophysical Journal Letters, 552(1):L1, 2001. 

[7] H. C. Chiang and et. al. Measurement of Cosmic Microwave Background Polarization Power 
Spectra from Two Years of BICEP Data. The Astrophysical Journal, 711(2):1123, 2010. 

[8] M. Cortes, A. R. Liddle, and P. Mukherjee. On what scale should inflationary observables be 
constrained? Phys. Rev. D, 75:083520, Apr 2007. 

[9] M. Cortes, A. R. Liddle, and D. Parkinson. On the prior dependence of constraints on the 
tensor-to-scalar ratio. Journal of Cosmology and Astroparticle Physics, 2011(09):027, 2011. 

[10] J. Dunkley and et. al. The Atacama Cosmology Telescope: Cosmological Parameters from the 
2008 Power Spectra. [arXiv: 1009.0866], 2010. 



- 13 - 



[11] F. Feroz and M. P. Hobson. Multimodal nested sampling: an efficient and robust alternative to 
Markov Chain Monte Carlo methods for astronomical data analyses. Monthly Notices of the 
Royal Astronomical Society, 384(2):449-463, 2008. 

[12] F. Feroz, M. P. Hobson, and M. Bridges. MultiNest: an efficient and robust Bayesian inference 
tool for cosmology and particle physics. Monthly Notices of the Royal Astronomical Society, 
398(4):1601-1614, 2009. 

[13] Z.-K. Guo, D. J. Schwarz, and Y.-Z. Zhang. Reconstruction of the primordial power spectrum 
from CMB data. Journal of Cosmology and Astro-particle Physics, 2011(08):031, 2011. 

[14] Z.-K. Guo and Y.-Z. Zhang. Uncorrelated estimates of the primordial power spectrum. Journal 
of Cosmology and Astroparticle Physics, 2011(11):032, 2011. 

[15] R. Keisler and et. al. A Measurement of the Damping Tail of the Cosmic Microwave 

Background Power Spectrum with the South Pole Telescope. The Astrophysical Journal, 
743(1):28, 2011. 

[16] E. Komatsu and et. al. Seven-year Wilkinson Microwave Anisotropy Probe (WMAP) 

Observations: Cosmological Interpretation. The Astrophysical Journal Supplement Series, 
192(2):18, 2011. 

[17] A. Lasenby and C. Doran. Conformal Models of de Sitter Space, Initial Conditions for Inflation 
and the CMB. AIP Conference Proceedings, 736(1) :53-70, 2004. 

[18] A. Lasenby and C. Doran. Closed universes, de Sitter space, and inflation. Phys. Rev. D, 
71:063502, Mar 2005. 

[19] A. Lewis and S. Bridle. Cosmological parameters from CMB and other data: A Monte Carlo 
approach. Physical Review D, 66(10), 2002. 

[20] A. Lewis, A. Challinor, and A. Lasenby. Efficient Computation of Cosmic Microwave 

Background Anisotropies in Closed Friedmann-Robcrtson- Walker Models. The Astrophysical 
Journal, 538(2) :473, 2000. 

[21] A. R. Liddle and D. H. Lyth. COBE, gravitational waves, inflation and extended inflation. 
Physics Letters B, 291(4):391 - 398, 1992. 

[22] Y.-Z. Ma, W. Zhao, and M. L. Brown. Constraints on standard and non-standard early 
universe models from CMB B -mode polarization. Journal of Cosmology and Astroparticle 
Physics, 2010(10) :007, 2010. 

[23] D. Parkinson and A. R. Liddle. Application of Bayesian model averaging to measurements of 
the primordial power spectrum. Phys. Rev. D, 82:103533, Nov 2010. 

[24] L. Perotto, J. Lesgourgues, S. Hannestad, H. Tu, and Y. Y. Y. Wong. Probing cosmological 
parameters with the CMB: forecasts from Monte Carlo simulations. Journal of Cosmology and 
Astroparticle Physics, 2006(10) :013, 2006. 

[25] Planck Collaboration. The Science Programme of Planck, [astro-ph/0604069]. 

[26] B.A.Powell. Tensor tilt from primordial b-modes. [arXiv: 1106.5059], 2011. 

[27] B. A. Powell. Scalar runnings and a test of slow roll from CMB distortions. [arXiv: 1209.2024], 
2012. 

[28] B. A. Reid and et. al. Cosmological constraints from the clustering of the Sloan Digital Sky 
Survey DR7 luminous red galaxies. Monthly Notices of the Royal Astronomical Society, 
404(l):60-85, 2010. 

[29] A. G. Riess and et. al. A Redetermination of the Hubble Constant with the Hubble Space 
Telescope from a Differential Distance Ladder. The Astrophysical Journal, 699(1):539, 2009. 

[30] R. Trotta. Bayes in the sky: Bayesian inference and model selection in cosmology. 



- 14 - 



Contemporary Physics, 49(2):71-104, 2008. 

[31] J. A. Vazquez, M. Bridges, M. Hobson, and A. Lasenby. Model selection applied to 

reconstruction of the Primordial Power Spectrum. Journal of Cosmology and Astroparticle 
Physics, 2012(06):06, 2012. [arXiv:1203.1252]. 

[32] J. A. Vazquez, A. N. Lasenby, M. Bridges, and M. P. Hobson. A Bayesian study of the 

primordial power spectrum from a novel closed universe model. Monthly Notices of the Royal 
Astronomical Society, 422(3):1948-1956, 2012. 



- 15 - 



